Magnetization switching in nanoscale ferromagnetic grains: 
simulations with heterogeneous nucleation 
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We present results obtained with various types of hetero- 
geneous nucleation in a kinetic Ising model of magnetization 
switching in single-domain ferromagnetic nanoparticles. We 
investigate the effect of the presence of the system boundary 
and make comparison with simulations on periodic lattices. 
We also study systems with bulk disorder and compare how 
two different types of disorder influence the switching behav- 



PACS Number(s): 75.50.Tt, 75.40.Mg, 64.60.Qb, 05.50.+q. 



This work is devoted to computer simulation of 
magnetization switching in single-domain ferromagnetic 
nanoparticles. For highly anisotropic systems, simple ki- 
netic Ising models cap, qualitatively explain many exper- 
imental observations a However, to obtain more realistic 
results additional physical effects must be included in 
the model. An important aspect of real samples is het- 
erogeneous droplet nucleation due to the presence of the 
system boundary and due to bulk disorder. 

The model used in our study is a square-lattice nearest- 
neighbor kinetic Ising ferromagnet with random up- 
dates using either Glauber or Metropolis single-spin-flip 
Monte Carlo dynamics. We first investigate systems 
with boundaries (but without disorder) , then we consider 
disordered systems with periodic boundary conditions. 
We concentrate on two quantities: the lifetime and the 
switching field. 

To determine the metastable lifetime, we start the sim- 
ulations with all spins +1 and impose a negative mag- 
netic field. Then we measure (in Monte Carlo Steps per 
Spin, MCSS) the mean time the sample needs to reach 
zero magnetization. Having measured this "lifetime" r 



for several values of the magnetic field, one can extract 
from the data the switching field H sw , which is defined 
as the field that produces a given lifetime. This quantity 
is often measured in experimentsBij 

We performed Monte Carlo simulation studies of the 
switching dynamics of a kinetic Ising model defined on 
circular lattices (which we define as subsets of a square 
lattice, contained within a circle of given diameter L). 
The Hamiltonian is 



H= - J OjOj - H 2J en - Hy, en 



(1) 



<ij> 



where the first two terms represent the standard spin- 
spin interaction with positive J and the coupling to the 
external field H, respectively. The last term, in which 
the summation runs only over the sites on the boundary 
of the lattice fi, is included to model the effect of the sys- 
tem boundary on nucleating droplets of the stable phase. 
For ifs = one has a free boundary, whereas positive 
(negative) He mimic a boundary which effectively repels 
(attracts) the droplets. 

Figure [l] shows the switching fields for a fixed wait- 
ing time r versus the diameter L of the circular lat- 
tice. For the smallest systems, H sw is zero; for slightly 
larger systems, H sw increases sharply with L. This is the 
coexistence region, in which the stable and metastable 
phases practically coexist.Era For larger systems, the be- 
havior depends strongly on the affinity of the bound- 
ary for droplets of the equilibrium phase, which here 
is modeled by different values of H^; and also on the 
waiting time r. In general, the increase in the coexis- 
tence region is followed by a decrease of the switching 
field in the single-droplet region (where switching is trig- 
gered by a single critical dropletcm). This results in a 
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maximum located at the crossover between the coexis- 
tence and single-droplet regions.cl Similar switching-field 
peaks are observed in certain experiments on nanoscale 
ferromagnets.0 For larger systems, more than one droplet 
nucleates during the switching process, and the system is 
in the multidroplet region EH There, the switching field 
becomes asymptotically independent of the system size. 
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FIG. 1. The switching field i/ sw vs. L for different 
values of the boundary field Hs for circular systems at 
T = 1.3J ps 0.57T C . The waiting time is r = 30000 MCSS, 
the dynamics is Metropolis. Data for periodic systems (full 
line) are shown for comparison. Note that even lattices with 
L — 400 are still in the single-droplet region. The decrease of 
the switching field would continue if we went to even larger 
systems. Eventually, for very large systems, it would converge 
to the same value in the multidroplet region, independent of 
He. 

Just as in periodic systems,! this behavior can be 
clearly observed in samples in which the boundaries re- 
pel droplets of the stable phase (positive Hs). For neu- 
tral boundaries (He — 0) the maximum in the switching 
field can be much less pronounced or it can disappear 
completely, depending on the waiting time. For exam- 
ple, with r = 1000, the maximum of H sw completely 
disappears for — 0. In general, to observe the maxi- 
mum, one needs a longer waiting time. Note that a long 
waiting time probably corresponds better to most real 
experimental situations. For theoretical considerations 
concerning magnetization switching in kinetic Ising sys- 
tems without disorder, see Refs. | and § 

To study effects of disorder, we simulated the above 
Ising model on periodic square lattices using Glauber dy- 
namics. The first type of disorder we investigate is pro- 
duced by defects generated by randomly deleting bonds 
of the lattice with concentration c. 

In order to understand how the disorder influences the 
switching, we measured various properties of the nucle- 
ating droplets of the stable phase during the simulation. 
All measured quantities were taken as mean values over 



all events in which a droplet of a given size appeared in 
the system. It turns out that with this type of disorder, 
the switching can be well described in terms of "average" 
nucleating droplets. 

Figure shows the critical size of a droplet versus the 
defect concentration c. Droplets larger than the criti- 
cal size are more likely to grow than to shrink, whereas 
smaller droplets do the opposite. As expected, the disor- 
der leads to smaller critical droplets, which in turn results 
in a dramatic decrease of the metastable lifetime. 
30 
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FIG. 2. Size and defect content of a critical droplet as 
functions of the defect concentration. 
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3. Switching field as a function of system size for 
different defect concentrations, waiting time r = 30000 MCSS 
and Glauber dynamics at T = 1.3J. 

Another interesting property of the typical critical 
droplet is its defect content: the number of deleted bonds 
associated with the droplet. It can be deduced from 
Fig. ||, as well as from direct measurements, that the 
effective concentration of defects within the droplets is 
considerably larger than the global concentration. This 
is caused by the fact that the droplets are preferentially 
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nucleated in the vicinity of the defects. 

In Fig. H we show the effect of the disorder on the 
switching field. Whereas the switching field is decreased 
by the disorder, its shape as a function of the system size 
remains approximately independent of c. 

Thus, in the intermediate-temperature region studied 
here, the switching dynamics of the kinetic Ising model 
on a bond-diluted lattice is essentially the same as in sys- 
tems without disorder, although the metastable lifetimes 
and switching fields are considerably reduced by the dis- 
order. The effect of the disorder can be approximately 
described in terms of an effective medium which affects 
the growing droplets ofrthe stable phase. Details will be 
given in a future paper B 
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FIG. 4. The metastable lifetime r as a function of dis- 
order. The error bars show the width of the distributions of 
lifetimes for different realizations of the disorder. For the di- 
luted-bond type of disorder (insert), this broadening is smaller 
than the symbol size. In both cases, T — 1.3J, H = — 0.3J 
and L = 20. 



of switching fields. The mechanism leading to this be- 
havior is that the size of a critical droplet is different in 
different parts of the disordered lattice. 

In summary, we have studied magnetization switching 
in systems in which heterogeneous nucleation of the equi- 
librium phase may occur on the system boundary or be 
associated with bulk impurities. 

The presence of the system boundary strongly affects 
the magnetization switching in small systems since the 
critical fluctuations of the stable phase tend to nucleate 
in its vicinity. This considerably decreases the metastable 
lifetime as well as the switching field, compared to peri- 
odic systems. However, the basic features, such as the 
crossover from the coexistence region to a single-droplet 
region, and then to a multidroplet region with increasing 
L, can still be observed in systems with a boundary. The 
switching field as a function of the system size becomes 
less peaked, but a maximum, located near the crossover 
between the coexistence and the single- droplet regions, 
can still be observed if the waiting time is sufficiently 
long. 

We have studied two different types of disordered sys- 
tems. Diluted-bond disorder can be a toy model for sam- 
ples with impurities, whereas the model with random 
magnetic moments could represent films with fluctuat- 
ing thickness. Our results show that whereas these types 
of bulk disorder result in shorter metastable lifetimes in 
general, the details of the magnetization switching can 
differ considerably for different types of disorder. 
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Another type of disordered kinetic Ising model is de- 
fined by the Hamiltonian Ti = —J'Yli<ij~ > a i a 3' J i (J j ~ 
H^2i a i a i, where the amplitudes Oj represent random 
local magnetic moments. Here, we take them uncorre- 
lated with a "box" distribution of width W, centered 
around unity. The switching dynamics with this type of 
disorder differs from the one discussed above in several 
respects. First, the approach in terms of "average" nucle- 
ating droplets is not useful because the interesting quan- 
tities have broad distributions. Second, there is a lack of 
self-averaging, even in the intermediate-temperature re- 
gion studied here. This is demonstrated in Fig. ^, where 
we show how the mean lifetime is reduced by the disorder. 
At the same time, the relative width of the probability 
distribution of the mean individual-sample lifetimes in- 
creases. Whereas switching of a sample with a particular 
realization of the disorder is more deterministic than in 
a pure system, an ensemble of disordered systems ex- 
hibits a wide distribution of lifetimes and, consequently, 
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